Modified brane cosmologies with induced gravity, arbitrary matter content and a 

Gauss-Bonnet term in the bulk 
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We extend the covariant analysis of the brane cosmological evolution in order to take into account, 
apart from a general matter content and an induced-gravity term on the brane, a Gauss-Bonnet 
term in the bulk. The gravitational effect of the bulk matter on the brane evolution can be described 
in terms of the total bulk mass as measured by a bulk observer at the location of the brane. This 
mass appears in the effective Friedmann equation through a term characterized as generalized dark 
radiation that induces mirage effects in the evolution. We discuss the normal and self-accelerating 
branches of the combined system. We also derive the Raychaudhuri equation that can be used in 
order to determine if the cosmological evolution is accelerating. 



I. INTRODUCTION 

The cosmology of the Randall-Sundrum (RS) model 
with a single positive-tension brane [1] is a viable pro- 
totype for the evolution of a Universe identified with a 
hypersurface in a higher-dimensional non-compact space. 
The crucial property of the model that leads to the emer- 
gence of a realistic low-energy evolution is the effective 
compactification of gravity around the brane within the 
AdS background. In the simpler version [2] the mat- 
ter is assumed to be localized on the brane, while the 
bulk space includes only a negative cosmological con- 
stant. The cosmological evolution can be identified with 
the motion of the brane within a static bulk [.]] . 

An obvious generalization of this picture takes into ac- 
count the possible presence of matter both in the bulk 
and on the brane, with the possibility of energy exchange 
[4, 5, 7]. It is remarkable that the general situation for 
an arbitrary bulk content has a very simple description 
[8, 9]. The effect of the bulk on the brane cosmologi- 
cal evolution, as determined by the effective Friedmann 
equation, can be incorporated in a single parameter char- 
acterizing the "strength" of the induced modifications: 
the integrated mass Ai of the bulk fluid as measured by 
a bulk observer at the location of the brane. If the spa- 
tial part of the brane has the geometry of a sphere, this 
mass is the same as the effective gravitational mass of 
the bulk. In the majority of cases the bulk observer is 
assumed to be comoving with the bulk fluid. For this 
reason we shall refer to the bulk mass as the comoving 
mass. We employ the same terminology for matter com- 
ponents, such as a bulk radiation fluid, for which there is 
no comoving observer. The comoving mass M in general 
depends on the brane scale factor R and the proper time 
T on the brane. It appears in the effective Friedmann 
equation within a contribution ~ Ai {R, r) / that has 
been termed generalized dark radiation [8, ')]. 

The next generalization includes terms that can be as- 
sumed to arise at the level of radiative corrections. The 
breaking of translational invariance by the location of 



the brane allows the presence of an induced gravity term 
on it. For a tensionless brane in a Minkowski bulk, one 
obtains the Dvali-Gabadadze-Porrati (DGP) model [10]. 
Specific examples of induced gravity can be obtained in 
string theory and are common in holographic descriptions 
[11, 12, 13] In the context of brane cosmology the pres- 
ence of the induced-gravity term has a remarkable con- 
sequence: the appearance of a self-accelerating branch in 
the brane evolution [14, 15]. The existence of matter in 
the bulk can be taken into account in complete analogy to 
the RS case [ I (>]. Its effects can be incorporated in the co- 
moving mass A4 of the bulk. Exotic modifications of the 
brane cosmological evolution may arise [IG]. Despite its 
very interesting properties, the self-accelerating branch 
is known to suffer from ghost-like instabilities [17, 18]. 

Radiative corrections in the bulk generate higher cur- 
vature terms. In particular, the Gauss-Bonnet (GB) 
combination is the leading bulk correction in the con- 
text of string theory [19]. The cosmological evolution 
in the presence of a bulk GB term (with or without in- 
duced gravity) has been discussed extensively [20, 21, 
22, 2.3, 24]. An interesting feature is that it is possible to 
embed the brane so that it has self-accelerating cosmo- 
logical expansion. Unfortunately, the branch of solutions 
that displays this behaviour is known to be unstable with 
respect to small perturbations [25] (see also [26, 27]). Fi- 
nally, the combination of both the induced-gravity and 
GB terms leads to a multitude of cosmological scenarios, 
even in the absence of matter in the bulk [2S]. 

We are interested in generalizing this framework in or- 
der to take into account an arbitrary matter content of 
the bulk. The effective action we consider has the form 



S = 




(1) 



In the first integral, —A is the bulk cosmological con- 
stant (in general we assume A > 0), Rabcd, Rab the 
Riemann and Ricci tensors, R the curvature scalar of the 
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5-dimcnsional bulk spacctimc with metric gAB, and 

^GB = — ARabR^^ + RabcdR'^^'^^ (2) 

the GB term with couphng constant a. In the second 
integral, V is the brane tension, g^p the induced 4- 
dimensional metric on the brane, g^ its determinant, i?4 
the corresponding curvature scalar, and fc the character- 
istic length scale of induced gravity. The pure gravity 
part of the action includes the standard Einstein term, 
along with terms that could arise through quantum cor- 
rections. The matter contributions are arbitrary, and 
the effective action incorporates possible quantum correc- 
tions in this sector. We assume, however, that the cor- 
responding energy-momentum tensor is consistent with 
the underlying geometry. As a result, the framework we 
consider is very general: It corresponds to a generic low- 
energy effective action in the Einstein frame, involving 
possibly a multitude of fields, and includes the leading 
quantum corrections. 

The Einstein Field Equations (EFE) take the form 



of a bulk cosmological fluid with a certain equation of 
state. The form of A4{R,t) is necessary for a detailed 
discussion of the cosmological evolution of the brane. On 
the other hand, the general properties of the evolution, 
such as the presence of acceleration, can be determined 
from the structure of the equations we shall derive. For 
this reason we shall not consider specific models in this 
paper. We postpone a detailed investigation of the role 
of the comoving mass in individual cases for future work. 
Furthermore and although one can relax the assumption 
of a Z2-symmetry [30] for the sake of simplicity, we main- 
tain the existence of the mirror symmetry around the 
location of the brane. 

Throughout this paper the following conventions are 
used: the pair {Ai , g) denotes the 5D bulk spacetime 
manifold endowed with a Lorentzian metric of signature 
(—,+,+,+,+), bulk 5D indices are denoted by capital 
latin letters A,B,... = 0,1, 2,..., 4, greek letters denote 
brane indices a, /?,... — 0,1,2,3, and lower case latin 
letters indicate spatial 3D components. 



aH\ = G\^^ (Tj, + AS^s) , (3) 



with the energy-momentum (EM) tensor given by 

TAB^Tfj^'' + S{v)TAB- (4) 

The corrections to the EFE (3) originating in the GB 
term are represented by the Lovelock tensor 

Hab = '2RRab ~ 4:Rj^^ Rkb — 4:R^^Rakbl + 



-2i? 



KLM 



RbKLM — -j^gAB^GB- (5) 



The term is the bulk matter contribution, while 

TAB is the contribution from the brane located at 
ri{x^) = and has the form [J:;] 



Tap 



T, 



Vg^p - 2r,Af 3 G 



Q/3 



(6) 



where T^^^^ is the brane Energy-Momentum (EM) ten- 
sor. We note that the presence of the induced 4- 
dimensional curvature term results in a contribution to 
the tensor Tap proportional to the Einstein tensor Gap 
on the brane. 

The purpose of this paper is to provide the general 
form of the solution of the above equations in the case 
of a FRW brane, for which there exist 3-dimensional hy- 
persurfaces T) invariant under a six-dimensional group of 
isometrics. It follows that the surfaces T) have constant 
curvature, parametrized by the constant fc = 0,±1. As 
we have already mentioned, the solution depends on the 
comoving mass Ai of the bulk fluid, which is a function 
of the brane scale factor R and the proper time r on the 
brane. The form of this function can be determined only 
within a specific model of the bulk dynamics. Such a 
model may involve several bulk fields that possibly inter- 
act with the brane, or may employ a description in terms 



II. 3-BRANE EMBEDDING IN A STATIC BULK 

Before presenting general and covariant results for the 
brane evolution, it is instructive to analyze a class of 
cases in which the problem is tractable in specific coor- 
dinate systems. We assume that for a certain observer 
the bulk content can be described as a static fluid. This 
assumption allows the possibility of an arbitrary number 
of fields and relies only on the existence of an observer 
comoving with the bulk matter. Clearly, important phys- 
ical situations, such as those that involve the propaga- 
tion of electromagnetic or gravitational radiation, are 
excluded by our assumption. However, many interest- 
ing backgrounds, including generalized black-hole ones, 
are allowed. 

In order for the embedding of a cosmological 3-brane 
to be possible, the spatial part of the metric must include 
a 3-space of constant curvature. The resulting metric can 
be cast in the form 



The Ihs of the EFE (3) take the form 
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(9) 
(10) 



where the prime denotes a derivative with respect to r. 



3 



The general form of the bulk energy-momentum tensor 
for the above geometric setup is 



with m(r, ?7 = 0) = 1. 
nate transformation 



Through an appropriate coordi- 



Tl'i3'^'' = diag(-p,p,p,p,p), 



(11) 



with the two pressures p, p not equal unless the bulk 
matter can be interpreted as a perfect fluid. The 00 com- 
ponent of (3) gives 



2a Akc 



1 



^-i + iFJ (12) 



whereas the combination of the 00 and 44 components 
results in 



{bn)' 
bn 



4a 



4ak\ _ 1 
~r) ^ GAP 



b\{p + p). (13) 



The conservation of the bulk energy-momentum tensor 
can be written in the form 

p' 3(p - p) _ (6A/3) -1 (j5 + xy^b'^ + krb^ ~ r 



p + p r{p+p) 



r2 - 4a&-2 4ctk 



(14) 

Because of the Bianchi identities, the set (12)-(14) com- 
pletely describes the solution. 
Integrating (12) we find 



2a 4ka 



kr 



M{r) 



62 b"^ b^ 12M^ Gtt^A'P 

where A4 (r) satisfies 



dM 

dr 



2tt'^Pp 



2afc2, (15) 



(16) 



and corresponds to the comoving mass of the bulk fluid. 
The above equation has the solutions 



1 

b^ 4a 



2a A 4a7W(r) 



3A/3 Stt^MV^' 



(17) 



with El ~ ±1. For a the solution with ei = 1 
reproduces the known expression in the absence of the 
GB term [5]. This expression includes a contribution 
~ Ar^ arising from the bulk cosmological constant. 

In the branch described by the solution with ei = —1, 
the leading contribution to 1 /b"^ for a — > is ~ r^/a. One 
expects behaviour similar to that arising from an effective 
bulk cosmological constant 1/a. A brane embedded in 
such a background can display self-accelerating cosmo- 
logical expansion with constant ^ 1/a. Despite its 
very interesting properties, this branch is known to be 
unstable with respect to small perturbations [25]. 

In order to analyze the cosmological evolution of the 
brane, we employ the Gaussian normal coordinate system 
in which the metric takes the form 



ds^ 



-m^iT, r])dT^ + a^{T, ifjdVil + drf 



(18) 



r{T,ri) 



(19) 



the metric (7) can be written in the form of equation 
(18). We define i?(r) = a(T,ri = 0). In the system of 
coordinates (t, r) of equation (7) the brane is moving, as 
it is located at r = R{t). Hence [-5] 



dt_ 
drj 
da 

da 
drj 



1 



n{R) 



b\R)R^ + 1 



1/2 



R 



b{R) L 



b'^{R)R^ + 1 



1/2 



(20) 
(21) 
(22) 
(23) 



where the dot denotes a derivative with respect to proper 
time and 62 = ±1. The 77-derivatives are evaluated for 
?/ = 0^. The value of £2 determines the way the brane 
is embedded in the bulk space. As we have mentioned 
earlier, we impose a Z2-symmetry around the brane. For 
a matter configuration that solves the EFE in an infinite 
bulk before the brane embedding, only the solution in 
the half-space and its mirror image are employed in the 
construction that includes the brane. The value of e2 
determines which half-space is used [18]. A negative sign 
in the r.h.s. of equation (23) means that r decreases 
away from the brane. In the absence of induced gravity 
and a GB term, the brane has positive tension. The 
configuration is stable under small perturbations and the 
massless graviton is localized near the brane. 

The bulk energy-momentum tensor at the location of 
the brane in the coordinate system (r, rj) is 



ry-iBULK 

^00 



rriBULK 

J 44 

thBULK 

-^04 



p{R) + [p{R)+p{R)]b'{R)R' 
R^p{R) (no summation) 

p{R) + [p(R) + p{R)]b\R)R^ 
1/2 



e2b{R)R ¥{R)R 



1 



(24) 
(25) 
(26) 

[p{R)+p{R)l27) 



The sign of ro4 '''^ indicates whether a brane observer de- 
tects inflow or outflow of energy. This sign is determined 
by the value of £2- For £2 = 1 the volume of the bulk 
space, as well as the matter it contains, grow for increas- 
ing rj. Conservation of energy requires that there is en- 
ergy outflow from the brane. For £2 = — 1 the brane 
embedding is such that the bulk volume diminishes for 
increasing 77. This is consistent with energy flowing into 
the brane from both sides. 

The Ihs of (3) near the brane [rj ^0^) take the form 
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(no summation over repeated indices) 



+- 



12aa"a'^ 3a^ 

a? 

3a" I2aka" 3fc 



3a2 



(28) 



4am" (a')^ [a'f 2m' a' iam'a"a' 



8a {m')^ Aam"d 
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2a" iahn" 

+ — 

a a'' 

&aa"d 2d 



m" 



16am'dd' 



(29) 



G\ 



12am' {a'y 12ad{a'y 3 (a') 



l\2 



+ - 



a-^ a"^ 
Uam'd'^a' Uakm'a' 3m'a' 3(2^ 



a 

3k 12ad^d 12akd 3d 
_ 



a 



a" 



'a2 



a 



(30) 



12am' 12ad'd2 12a(a') m'd 
1 1 — 

a^ a^ a^ 

\2akm'd 3m'd 12a (a')^ d' 



+- 



12afcd' 3d' 



(31) 



with a prime now denoting a derivative with respect to 
77. 

We consider a brane Universe containing a perfect fluid 
with an energy-momentum tensor 



rriBRANE X/'„\„2 



= 5{ri)a (r, 77)diag 



a^(r, rf) 



P,P,P,P,0 



(32) 



Integrating the 00 and ii components of (3) on a smah 7/ 
interval around the brane and using (6) we obtain 



1 + 4a 



12M3 




V + p- 6rcM^ 



(33) 



m 



1 + AaH-" 



Aka 

r,2 



Aaa' 



I 2 



4M3 



P-V + 2rcM^ H 



2d 
a 



(34) 



From (33) and (23) it is straightforward to derive the 
effective Friedmann equation 



1 



1 



144M6 



b'^a^ 
V + p 



1 



1 2 



In the low-energy limit (p, iJ, a 



(35) 

0) our choice of sign 



for a'j^ in (23) must be consistent with the rhs of (33). 
For example, for a = Tc = a negative sign (£2 = 1) 
for a'_^_ is consistent with the negative sign in the rhs of 
(33) only if F > 0. It is worth mentioning that the 
above choice is referred as the "normal branch" and con- 
tains the Randall-Sundrum model as a particular case. 
In addition, the positivity of the brane tension guarantees 
stability under small perturbations and graviton localiza- 
tion near the brane. The self-accclcrating branch of the 
DGP model [10, 14, 15] has a = V ^ 0, p,a-^ ^ and 
^ l/vc- It is then apparent from (33) that a(|_ > 
(e2 = ^1)- This branch is known to have ghost-like in- 
stabilities [17]. The effective Friedmann equation (35) 
results from squaring a'_^_. As a result, it includes both 
the normal and self-accelerating branches of the DGP 
model. 

We define the effective cosmological constant 



A = 



4M6 



1 - Vl - A 



^ 2 



(2 + v^l - a) 



(36) 



where A = 2aA/3M3. The low-energy cosmological evo- 
lution can be determined by expanding the Friedmann 
equation (35) in a~^ and p. We assume that the cosmo- 
logical constant has been tuned to zero (A = 0). In the 
normal branch (no self-acceleration) with ei = e2 = 1 we 
expect the standard behaviour ~ yO for y 7^ 0. The 
Friedmann equation becomes 



+ — ^ (72M^ + IGaAM^ + GrcVM^Y 



-,(2+v/r^)M3-^ 



5 7W(a) 



Vp 



(P + Prf). (37) 



The effective Planck constant is 



and the generalized dark radiation [9] 



M^i = ( 1 



2aA r,V \ 12M^ 

9AP + 12A/3 ) V ^ ^ 



Pd 



2+ v/l-A12M3A4(a) 



(39) 
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The latter quantity desribes a mirage energy density that 
affects the evolution without arising from a source on the 
brane [31]. In the last two expressions we have absorbed 
the dependence on a and Vc in terms that approach 1 
for Tc, a — *■ 0. It is remarkable that, even in the pres- 
ence of GB and induced-gravity terms, the effect of the 
bulk matter in the low-energy limit can be absorbed in 
a mirage density term ^ A4{R) / R^. More exotic low- 
energy behaviour can be observed in the self-accelerating 
branch. 

A simple example of the cosmological evolution we de- 
scribed can be obtained if the bulk energy-momentum 
tensor is assumed to have the form (11) with p = p. In 
order to avoid the appearance of metric singularities in 
equation (15) we assume that k ^ 1. If the equation of 
state p = p{p) is known, the bulk metric can be com- 
pletely determined. This background is a generalization 
in four spatial dimensions and for a negative cosmological 
constant of the conventional solution describing the inte- 
rior of stars. For this reason it has been termed AdS-star 
in reference [■">]. 

The presence of a GB term leads to quantitative mod- 
ifications of the a = configuration, but its qualitative 
form remains the same. Moreover, for an equation of 
state of the form p = wp"^ the asymptotic form of Ai{r) 
for large r is the same for all values of a. This means that 
the late-time cosmological evolution is given by equations 
(37)-(39) with M{R) as in reference [5]. 



III. COVARIANT STRUCTURE OF THE 5D GB 
BULK WITH A FRW BRANE 

As was demonstrated in [8] the bulk geometry with a 
FRW brane can be seen as a 5-dimensional generalization 
of the inhomogeneous orthogonal family of Locally Rota- 
tionally Symmetric (LRS class II) spacctimes [.'J]. The 
orbits V of the six-dimensional multiply transitive group 
of isometrics are maximally symmetric 3-dimensional hy- 
persurfaces with spatial (constant) curvature determined 
by the value of fc = 0, ±1. The rotational symmetry of 
the bulk manifold Ai results in several key features of its 
geometric structure that allow a simplified and unified 
treatment. To begin with, let us first note that, from a 
dynamical point of view, there arc two different ways to 
choose the unit timelike vector field normal to the group 
trajectories: either adapted to the average fluid velocity 
of the bulk matter configuration or to the prolongated 
brane observers u^. In coordinate language, this freedom 
is related to the choice of a particular coordinate system. 
As we are interested in a bulk that is not empty, it is 
convenient to use the velocity of the bulk observers 
{u-^ua = ~1) in what follows. 

Using the standard 1-1-4 splitting of the 5D spacetimc 
manifold the deformation of the timclikc congruence 
can be expressed in terms of the corresponding kinemat- 



ical quantities 

UA-B = ^AB + ^^AB ~ UAUB = VAB - UAUb, (40) 

where Y^ab = {^^^g - \h^^hAB) U(k-l), © = w^, and 
ii'^ = u^u^ are the rate of shear tensor, the rate of 
expansion scalar, the vorticity tensor and the acceleration 
of the observers m"*, respectively, Hab = Qab + uaub is 
the projection operator perpendicularly to u^, and vab 
is the extrinsic curvature of the 4D spaces S normal to 
u^. We recall that the fluid velocity is orthogonal to 
the group orbits, so that the vorticity tensor is ^ab — 0. 
This fact allows us to interpret Hab as the metric of S 
and employ an appropriate covariant derivative inherent 
to these spaces 

D,P^^ -,j,„ ^ hih^s-h]hj...hf (P^' -rx...).^^ (41) 

for any tensor P^^ - jj . The structural characteristics 
of iS are described in terms of the kinematical quantities 
of u^ by using the Gauss equation 

^RaBCD = h/hB^hc^'hj^RKLMN + 2vA[DVc]B, (42) 

with *Rabcd the curvature tensor of S. 

The assumption of maximal symmetry of the group or- 
bits V implies the existence of a preferred spacelike direc- 
tion {e^BA = 1, u^ca = 0), that represents the local 
axis of symmetry with respect to which all the geometri- 
cal, kinematical and dynamical quantities are invariant. 
As a result, all the spacelike vector or traceless tensor 
fields which are covariantly constructed via the timelike 
vector field u"^ can be written in terms of [n]. In or- 
der to study the structure of the spacelike congruence of 
curves generated by the unit spacelike vector field we 
proceed in complete analogy with the 1-1-4 decomposi- 
tion. The starting point is to introduce the projection 
tensor: 

IIab = 9AB + uaUb - cacb = hAB - bacb (43) 

n/ = 3, Uc^Iljf = n/, n/e^ = U/u^ O (44) 

which is identified with the associated metric of the 
3D manifold V (the screen space) normal to the pair 
{^^,6'^} at any spacetime event. The geometric struc- 
ture of V is analyzed by decomposing into irreducible 
kinematical parts the first covariant derivatives of the 
spacelike vector field according to [•■).■]] 

eA;B = Tab + -^^AB +TIab + SaCb ~ caub + 

+nB^ecUA + [2ncBe^ - Nb] ua- (45) 

Here 

d = eA;Bn^^ = e^.A + e^^A (46) 
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Tab - n/n 



Ta-lH^^ = 



(47) 
(48) 

(49) 



are the rate of the surface expansion, the rate of shear 
tensor, the rotation tensor and the Greenberg vector field 
of the spacehke congruence e"^, respectively. We use the 
notation 



K' =KA..,Le 



(50) 



for the directional derivative along the vector field of 
any scalar or tcnsorial quantity. 

Each of the above kincmatical quantities carries infor- 
mation on the (overall or in different directions) distor- 
tion of D as measured by the bulk observers u^. They 
have a similar interpretation as the corresponding quan- 
tities of the timclikc congruence u^. The new ingredient 
is the Greenberg vector Na which represents the "cou- 
pling" mechanism between directions normal and paral- 
lel to the screen space T>. For example, the equation 
N"^ = implies that the pair of vector fields {u^,eA^ 
generates a 2-dimensional integrable submanifold of M 
and the spacelike congruence is "comoving" ( "frozen- 
in") along the worldlines of the fundamental observers 
u'^. In addition, it ensures that Tab and TZab lie in 
the screen space and the unit vector fields |e'^,w'^} arc 
orthogonal at any instant. 

An important consequence of the preferred spacelike 
direction, or equivalently the induced three dimensional 
isotropy, is the fact that any, covariantly defined via u^, 
spacelike and traceless tensor field lying in the screen 
space must vanish. This means that = = Hg'^ec 
and Tab = = TZab-, and that the first derivatives of 
the spacclike congruence take the form 

eA-B = 3^^^ + e'AeB - caub = Kab + s'a^^b, (51) 



n4^n^'^e(/.j) is the extrinsic curvature of 



where Kab — ^^a 
the spacelike hypersurfaces S 

The vanishings of the Greenberg vector and the vor- 
ticity tensor imply that T) is an assemply of 3D hyper- 
surfaces that mesh together to generate the integrable 
manifold 2?, which is a submanifold of the observers' in- 
stantaneous rest space. Consequently, the fully projected 
(perpendicular to the pair {u^, e^}) covariant derivative 
"II" defined as 



P 



AB... 



IJ...\\L 



n/n/ 



•n/n/. 



RS.. 



TX.. 



■K ' 

(52) 

represents the proper 3D covariant derivative, since 
nAB||c = and A||[al] — foi' any scalar quantity A. 

The definition of the overall expansion d of the space- 
like congruence permits us to introduce the quantity £ 



according to 



J' 



\\A 



(53) 



Equation (53) makes clear the geometrical role of £ as 
the average length scale of V. For example, in the spher- 
ically symmetric case fe = 1 it represents the radius of the 
spheres V. On the other hand the temporal (it— )change 
of £ is controlled by the expansion rate of the timelike 
congruence as measured in the screen space I?, namely 



(54) 



Taking into account the above considerations, the length 
£ completely determines the volume of V which scales 
~ as the screen space V evolves. 

Regarding the dynamics, the matter content of the 
bulk is described by the energy- momentum T^j^^, which 
can be written in the usual way with respect to the ob- 



servers u 



Tab'^ = PUAUB + phAB + 2q(AUB) + TTAB- (55) 

The dynamical quantities measured by the bulk observers 
arc defined as 



rpHll-LK A B _ 1^ 



BULK 7,j4B 



AB 



, QA 



7 CtiBULK D 
'^A^CD " ' 



TTAB = h^Ah'^T^T - -{h^^T^''^^)hAB. (56) 

Because of the specific geometrical background of the 
bulk, it will be helpful to investigate the influence of 
the matter content on the curvature of the 3-dimensional 
screen space V. This can be achieved by using the fact 
that the screen space V forms an integrable submanifold 
of M with a well defined metric Hab and a proper co- 
variant derivative "||". Then the corresponding Gauss 
equation for the distribution normal to the 1— form u 
reads 

^RaBCD — ^A^B^C^^d' '^Rl.lKL + '2Ka[cRd]b, 

(57) 

where ^Rabcd is the curvature tensor of the screen space 
T>. Contracting twice equation (57) and using (42) we 
obtain 



_ 1 
T ^ 6 



H^^'U^^RabCD { ^Tl^'^UA-.B 



1 



(58) 

or equivalently in terms of the average scale factor and 
the 5D Einstein tensor, 

i^4(^ + iG.«."-2G.)- 



+ 



£' 



(59) 
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Here £ = C ac b dw^ is the spatial eigenvalue of 
the electric part of the Weyl tensor and SGj, = H^^Gab- 
Equation (59) shows how the scalar curvature of the 3- 
space T) is affected by the kinematics and the dynamical 
(when the 5D EFE are employed) content of the space- 
time. Equivalcntly it represents the evolution equation 
of the average length scale. 

We point out that equation (59) is a first integral of 
the propagation equation (along u^) of i/d, or the spatial 
variation (along e^) of the spacelike expansion In 
the context of brane cosmology, the physically interesting 
quantity is i/ i which {on the brane) corresponds to the 
overall expansion H of the 3-brane. The expansion of the 
timelikc congruence is written 



e 



1. 



-T^Ase e 



A„B 



(60) 



Then, the temporal projection of the trace and traceless 
symmetric part of the Ricci identities for gives evolu- 
tion equations for O and Y,ab- Combining the resulting 
expressions with the propagation of (60) we get 




-jeAU^ 



-Gabs^g 



+ 



\U + \gab9^''-2G, 



.(61) 



Obviously, equation (61) is equivalent to the Raychaud- 
huri equation for the expansion H of the brane Universe. 

If the GB correction is absent [a = 0) the first term in 
the rhs of equation (59) can be determined by using the 
full 5D EFE (3) and the bulk energy momentum tensor 
(55). The rcsuh is [8] 



M 



A 




6M37r2^4 12M3 



-CAU 



(62) 



M 



A 



6M3 



P\\ 



(63) 



where p\\ = T^'^^^e e is the pressure in the direction of 
the preferred spacelike vector field and A4 the comoving 
mass of the bulk fluid, satisfying 



{M - Mo)' = 2T:^pi^£'. 



(64) 



Only in the spherically symmetric case (fc = 1) the co- 
moving mass M has the usual physical interpretation of 
the effective gravitational mass contained within a spher- 
ical shell with radii £o and i. However, we shall refer to 
A4 as the integrated mass for all geometries of the hyper- 
surfaces 2?. The integration "constant" A4o in equation 
(64) can be interpreted as the mass of a black hole at 
4 = 0. 



When a 7^ the methodology breaks down because 
the EFE (3) include the Lovelock tensor as an additional 
contribution. Nevertheless, one can view 



Tab — 



1 



2A/3 



rriBULK Tj 

^ AB —aHAB 



(65) 



as an effective "energy-momentum" tensor, so that the 
EFE take their standard form and the methodology of 
[ ■] can be applied. This leads to 



1 



1 



2M3 



T-4pj 



where 



M 



2M^ai^^CGB) 



p^=p^- 2M''a-Il^^HAB 



(66) 



(67) 



(68) 



(^M ^ Mo)' = 27T^ {Ta 



A„,BijZiji 



B - 2M''aHAB) u'^u 

(69) 

It follows that the induced evolution of the average length 
scale i is 



k 



M 



H 



A 



M 



3^4 



where 



12M3 6Af37r2£4' 
(70) 



Mh = 2Habu 



'-HABU^U^i'^^ 



(71) 



Even though we have managed to express the eigenvalue 
of the electric part of the 5D Weyl tensor explicitly in 
terms of the "dynamical" variables of the effective bulk 
energy-momentum tensor (65), the appearance of the 
Lovelock correction term Mh makes the interpretation 
and the analysis of (70) unclear. It should be noted that, 
even though A^^f is defined in a similar manner as the 
comoving mass M of the bulk fluid, the character of the 
former is purely geometrical since it contains only terms 
quadratic in the 5D curvature. In particular, because of 
equations (42) and (57) the 5D Rabcd, Rab and R can 
be expressed in terms of the 3D curvature quantities, the 
extrinsic curvatures vab , Kab and their first derivatives 
along e^. Consequently, the Lovelock tensor is written 
as a quadratic expression involving '^R, vab, Kab- (We 
recall that the screen space V is maximally symmetric, so 
that the scalar '^R completely determines the curvature 
tensor.) 

After a tedious calculation and using equation (71) we 
get 



M 



H 



k 

£2 



(72) 
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Defining the quantity 



It follows that 



A 



we finally obtain 



k 

£2' 



(73) 




A 



6A/3 



aH-^ {e'^A'^y . (78) 



A + 2aA^ = - 



M 



12A/3 6Af37r2£4- 



(74) 



Clearly, the last equation coincides {off the brane) with 
the equation of motion of the sheets of the bulk matter 
configuration in the 5D GB gravity. The solutions of the 
quadratic equation (74) are 



A 



1 

4a 



4a 



2aA 



SAP 3A/37r2£4 



1/2 



(75) 



where ei = ±1 is the same as the one defined in equation 
(17). 

We conclude this section by noticing that with the 
above identifications, equation (63) is written as 




-^CAU^ -A- 



1 



GAP 



Similarly as before, we can argue that the quantity 
Habs^g^ has a geometric nature and should be express- 
ible in terms of '^R, vab, Kab- After a lengthy calcula- 
tion, we find a remarkable relation connecting this quan- 
tity with the evolution of the comoving mass, namely 



2HABe^e^£^H ■ 



-6 {i^A^y 



(77) 



IV. THE EFFECT OF THE BULK ON THE 
BRANE EVOLUTION 

In order to derive the brane cosmological evolution one 
must take into account the Israel-Darmois [34] junction 
conditions for the extrinsic curvature of the brane. As 
we have already mentioned, we impose a Z2-symmetry 
around the location of the brane. We assume that the 
brane Universe is filled with a perfect fluid, so that the 
brane energy-momentum tensor takes the form 



T, 



(79) 

where h^fj = gap+UaUfj is the projection tensor normally 
to the brane velocity Ua- We note that on the brane the 
metric of the 3D screen space V coincides with ha0 so 
that liap ~ hap and H = l/l represents the Hubble 
parameter. 

Essentially, equation (73) corresponds to the general- 
ized Friedmann on the brane. From this equation, it can 
be seen that the discontinuity enters as the first derivative 
of the average length scale along e'^ and is represented 
covariantly by the expansion of the spacelike congruence 
■d = Il"^ea-i3 = 3i'/£. Therefore, one should consider the 
junction conditions that involve only the spatial expan- 
sion i.e. the fully 11— projected of the first derivatives of 
e^. 

The covariant form of the junction conditions for 
braneworld models with a GB term in the bulk has been 
derived in [22]: 



Kaf3 + -^[9Jq/3 — 2Jgal3 — 2 {SPajpS + gapG-ys) K'^^] = ~ 



Tap - -^rgap 



(80) 



where 



Jap 



h2KKa^K\ + K^sK^^Ka^ 



2Ka^K''^Ksp - K^K, 



ap) 



(81) 



and (73), we get 



PaPyS — Rap-/5 + 2ga[s'^ R-f]P + 2gp['y'^ R5]a+ *Rga['yg5]p- 

(82) 

Bearing in mind that Gapu^'u^ = 3{H^ + k/i'^) [IG], 
contracting (80) with IIq,^, and using equations (6), (51) 




k 



12A/3 



V + p- 6r,AP [H +-p 



(83) 
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or, equivalcntly, 



£2 



1 + 



k Aq 
£2 



1 



12M3 



(p + V)- GrcAr H 



+ 


k 






-)] 




P ) 





1/2 



(84) 



where £2 = il depends on whether the spacehke expan- 
sion (equivalently ^' / (■) is negative or positive in the 
vicinity of the brane. This is the same parameter that 
appears in equation (23). The quantity ^0 = -^{tt^) is 
given by (75) and incorporates the effects of the bulk fluid 
[9]. Equation (84) is the generalization of the Friedmann 
equation on the brane given in [23]. As expected, the 
black hole mass A^o has been replaced by the comoving 
mass M. of the bulk fluid [')]. 

Although equation (84) is a third-order polynomial in 
B ~ + k/l"^, only one root is free of instabilities 
and reduces to the standard Randall-Sundrum solution. 
This can be seen by considering the pure induced-gravity 
model for which a = {). It follows from (84) that [f G] 



2 



k 

£2 



-£2 



Tc {V + ~P) 



Tc {V + P) 

12AP 
_^ rlM 



= 1 



^A 



1/2 



(85) 



The branch with £2 = — 1 reduces to the self-accelerating 
branch of the DGP cosmology [fO, 14, 15] for F = 0. 
This branch is known to suffer from instabilities under 
small perturbations [17]. The value £2 = 1 reproduces 
the normal stable branch of brane cosmology. For a 7^ 0, 
there are two possible values of Aq, given by equation 
(75). The solution with £1 = — 1 reproduces the self- 
accelerating brane cosmology in the presence of a GB 
term. However, the bulk configuration is unstable in this 
case [25]. 

For completeness, we also give the Raychaudhuri equa- 
tion in the presence of induced-gravity and GB terms. It 
follows from equation (78) and reads 



A 



6M3 



1 



6A/3 



where 



q = l/i={i/i)' 



(86) 



(87) 



denotes the acceleration parameter. We note that, for 
arbitrary bulk matter configurations, the comoving mass 
depends on both the average length scale and the time 
scale defined by the brane observers u^. However, for on 
brane considerations, we have ^o(^) and the acceleration 
parameter is written 



A 
6M3 



-Ao' 



1 

6M3 



pr 



d {tA, 
df 



(88) 



The junction condition for the discontinuous quantity 
cau^ follows from (80) 



i^-q^AB-{eAU^)+A(-^] {CAU^^ 



+4 - ('^^"^) -^^\P-^^ {2q + B)] . 

(89) 



We set 



A 



dl 



(90) 



where p|| = T^^^^e^e^ is the pressure in the preferred 
spacelike direction. By using equations (73), (83) and 
(86) we find 

(91) 

V2 

where 

Qi = 128ylga(2ra-l) + 32ylo[8ra(26a+l)- 

-ABa - 3] + 16r (86a + 3) + 46^ (64a - Sr^) + 
4y - 3p + 



+2B I 48 -h r, 
1 



g2 = 4 



{~p + V)(j>^V) 

MAla^ + ■i2Aoa - 256B^a^ + 
P + V 



-2B 3ri 



64a - 12 - 



Af3 



(92) 



(93) 



Starting from equation (91) we can check which values 
of the parameters of the theory can lead to accelerating 
expansion. For example, if we restrict our considerations 
to the normal branch with £1 = £2 = 1, both the induced 
gravity model (a = 0, Vc ^ 0) and the GB model (a 7^ 0, 
Tc = 0) require the standard mechanisms of inducing ac- 
celeration: either an effective cosmological constant, or 
a negative integrated mass of the bulk fluid, or a fluid 
with sufficiently negative pressure [9]. The full analysis 
of the various cosmological flows for all possible values 
of a, Tc, £1, £2 is beyond the scope of this work. For an 
AdS-Schwarzschild bulk, various possibilities have been 
studied in the literature. (See for example [28] and ref- 
erences therein.) 

One significant property of the solutions with a non- 
zero bulk integrated mass M is that the brane energy 
density has extrema. This can be seen by considering p 
as a function of B and ^0 



p + V = 2M^ 3rcB + 2e2\/ B - AoiSaB+AaAo + 3) 



(94) 
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We emphasize that if = the value of is fixed by 
A and a (see eq. (75)), so that it docs not vary with 
the scale £ and the following analysis is not valid. There 
exist extremal values 



obtained for 

B 



APtc (-2r2e2 + 48a - Sr^) 



— 16a 
64a2 ' 



32a2 



^0 



4a 



, (95) 



(96) 



For £2 = 1 we get a global minimal value for the brane 
density 



APrc (5r2 - 48a) 



(97) 



For 62 = —1 we get a global maximal value 

M^rc (rl - 4:8a) , , 

(P + nnax = ^- (98) 

If we require a positive brane tension, it is clear that 
we must impose r\ — 48a > for a physically relevant 
solution. 



V. CONCLUSIONS 

The three equations (75), (84), (91) are the main re- 
sults of our work. They determine the brane evolution 
for a given, but otherwise general, distribution of bulk 
matter. At the practical level, they demonstrate how to 
construct cosmological brane models in non-trivial bulk 
backgrounds. One first has to find a solution of the EFE 
in the bulk. The embedding of the brane is then auto- 
matic and leads to a cosmological evolution described by 
(75), (84) and (91). It must be emphasized that, from 
the point of view of the brane observer, the procedure de- 
scribed in the present paper leads to a definite prediction 
for the energy exchange rate. Assigning a physical mean- 
ing to this rate may not be always straightforward. In 
realistic cases the bulk configuration must have a certain 
freedom of parameters in order to accomodate the ex- 
change rate that is expected on physical grounds. A typ- 
ical example is provided by a brane that radiates or ab- 
sorbs massless particles (Kaluza- Klein gravitons or gauge 
fields). The bulk metric is of the generalized Vaidya type 
and the energy-momentum tensor contains a radiation 
field. This metric includes an arbitrary function that can 
be fixed by requiring a rate of energy exchange consistent 
with the physical processes assumed on the brane (e.g. 
energy collisions in a hot plasma) [(>]. We postpone the 
discussion of these issues in specific models for a future 
publication. 

The possibility of accelerated expansion of the brane 
Universe is related to the structure of the Raychaudhuri 
equation (91). Its complicated nature does not permit 



a straightforward analysis. However, one can still draw 
some intuitive conclusions regarding the initial state of 
the brane Universe. Using equation (94) the expression 
for the acceleration parameter takes the simple form 



q = {4M3 [2A (2ra - 1) + F + 2S] + 
+62 VS- A {2r,M'^B + p-V)}x 



4M3 ( AAaa - 8Ba - 1 - rce2 VB^^i 



(99) 

For the values of B and that correspond to the 
extrema of the brane energy density (equations (96)) the 
acceleration parameter becomes 

32M^rca + [APrc (r^ - 16a) - S2a^ {V - p)] £2 

^ ~ 128APrca^{l + 62) ' 

(100) 

In the branch with £2 = ~1 the brane energy density 
has a maximal value given by equation (98). The stan- 
dard picture of an initial singularity, accompanied by 
infinite energy density and deceleration, is replaced by a 
state of maximal energy density. However, from equation 
(100) we deduce that the acceleration parameter is infi- 
nite for any brane equation of state apart from p = —p. 
This signals the presence of a curvature singularity (see 
e.g. [28] for the simple case of a Minkowski bulk). On 
the other hand, if initially p = —p, using (95) we get 



16a 



64a2 



(101) 



For a > and under the constraint — 48a > imposed 
at the end of the previous section, the brane evolution 
begins with finite energy density and pressure, as well as a 
positive and finite acceleration parameter. We point out 
that the elimination of the initial curvature singularity in 
the £2 = — 1 branch depends crucially on the assumption 
Ao 0. This is only an example of the plethora of new 
phenomena that emerge from the combination of induced 
gravity on the brane, and a GB term and matter in the 
bulk. 
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